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Abstract 

We re-examine the experimental results for the magnetic response function x "{h E,T), for q around 
the anti-ferromagnetic vectors Q, in the quantum-critical region, obtained by inelastic neutron scattering, 
on an Fe-based superconductor, and on a heavy Fermion compound. The motivation is to compare the 
results with a recent theory, which shows that the fluctuations in a generic anti-ferromagnetic model for 
itinerant fermions map to those in the universality class of the dissipative quantum-XY model. The quantum- 
critical fluctuations in this model, in a range of parameters, are given by the correlations of spatial and of 
temporal topological defects. The theory predicts a x"(q, E, T) (i) which is a separable function of (q — Q) 
and of {E,T), (ii) at crticality, the energy dependent part is oc tanh(if/2T) below a cut-off energy, (hi) the 
correlation time departs from its infinite value at criticality on the disordered side by an essential singularity, 
and (iv) the correlation length depends logarithmically on the correlation time, so that the dynamical critical 
exponent z is oo . The limited existing experimental results are found to be consistent with the first two 
unusual predictions from which the linear dependence of the resistivity on T and the T In T dependence of 
the entropy also follow. More experiments are suggested, especially to test the theory of variations on the 
correlation time and length on the departure from criticality. 
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I. INTRODUCTION 


Soon after the discovery that in the normal metallic region for dopings near those for the 
high superconducting transitions in cuprates have transport and thermodynamic properties [I] 
unlike those expected of a Fermi-liquid, it was discovered that the heavy-Fermion compounds, 
for compositions near where their anti-ferromagnetic (AFM) transition temperatures T/v 0, also 
have similar anomalies [2]. For example, in both cases the resistivity has a temperature dependence 
proportional to T at low temperatures and the entropy or thermopower (entropy per carrier) has 
a contribution proportional to TlnT. More recently, the AFM quantum-critical region in the Fe 
based superconductors also show the same anomalies in the resistivity and the thermo-power (see, 
e.g.. Refs 0 and m). 

The theoretical study of quantum-criticality in AFM’s began long before these experiments 
came along, starting with the work of Moriya, and Hertz and others [5H8]. These are extensions of 
the theory of classical dynamical critical phenomena [9] to the quantum-regime and may be called 
renormalized spin-wave theories. Since the above experiments, an enormous amount of theoretical 
work in further developing theories based on similar ideas has been done [7] . The results from such 
theories appear to correspond to experiments in AFM’s in which the fluctuations are 3-dimensional 
(3d) |in( ITT] but not for those in which they are two-dimensional (2d). In the past few years, S- 
S. Lee and collaborators [12] have shown that such theories are not controlled in 2d. This has 
been followed by significant theoretical work which seeks to approach the actual 2d problem by 
expanding about limits HMD in which the theory is controlled, for example about a Fermi-surface 
in 3d while the spin-fluctuations are 2d. Another set of ideas for the heavy-fermion problem relies 
on the fact that an isolated Kondo impurity in a metal has singular properties near the criticality 
of the host m and by the approximation that the heavy-fermion metal may be regarded as a 
self-consistent set of periodic Kondo impurities using methods such as the Dynamical mean-Field 
Theory. Yet another seeks to understand the results within the renormalized spin-wave framework 
by invoking a phenomenological renormalization of the spatial correlation length m- More exotic 
ideas have also been proposed m- 

A radically different solution to the 2d AFM problem has recently been proposed [T8|. It 
starts by showing that the model for criticality of itinerant AFM’s is the dissipative XY model 
(with appropriate lattice anisotropy). This is true for the AFM with planar order, either about 
a commensurate or incommensurate wave vector, or Ising order about an incommensurate wave 
vector. The quantum-dissipative XY model [19] in 2d in such a region has been transformed [20] 


2 


to a model in which the critical properties (on the disordered side) are determined by topological 
excitations - 2d vortices as in the classical 2d-XY model and a new variety of topological excitations, 
called warps which are instantons of bound monopole-anti-monopoles with net charge 0. Two 
dimensionless parameters specify a line along which the quantum-critical point occurs in the model, 
which may be taken to be \/KKr and a. K is the Josephson coupling energy of the phase variables, 
K~^ is the magnitude of the kinetic energy parameter, and a is a dimensionless parameter for the 
dissipation due to decay of phase fluctuations [HI |22]. The quantum-criticality in this space, from 
the disordered state to the ordered state, appears in quantum Monte-Carlo calculations to have 
either a dynamical critical exponent z = 1 or z = oo, depending on the critical value of the 
parameters [23]. Only the z = oo case, which occupies the bulk of the parameter space in the 
— a plane is summarized below, as this alone provides results in correspondence which bear 
relation to the theory. 

A stringent test of the theory requires measurements of the absolute magnitude of x" 
in the fluctuation regime over a range of q and spanning frequencies from well above to well below 
T. The thermodynamic and transport properties can usually be derived from x"{^.^E,T). The 
purpose of this paper is to present the limited available existing quantitative experimental results 
for this function, which are available in the range of q, E, and T to compare and check whether 
there is a correspondence between the predictions of the theory. 

II. CORRELATION FUNCTIONS 

The correlations of the dissipative quantum 2d-XY model have been recently investigated in 
analytical calculations [20] . These have been checked and extended by quantum Monte-Carlo 
calculations [22] . These calculations confirm that the transitions are driven by vortex and/or warp 
binding and not by (renormalized) spin-waves which only serve to generate effective interactions 
among the topological charges. The calculated correlation functions as a function of distance r 
and imaginary time r on the disordered side by 

X(r,r) xo-e“^/^ln (—(1) 
r r' 

r Pc 

ir = \tc\ eV ^ ln(^^/|rc|). (2) 

Here r is the imaginary time, periodic in l/{kBT) [ks = 1 in this paper), which has a lower cut-off 
\tc\ ~ and rc is the short distance spatial cut-off, of the order of a lattice spacing, p is 

the parameter which tunes the transition by tuning KK^ or a, and Pc is the critical value of p at 
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which the transition occurs at T —)• 0. xo together with the high frequency cut-off serves to give 
the integrated magnetic fluctuations through the sum-rule on them. 

To compare with experiments, it is necessary to Fourier transform the correlation function to 
momentum (q) and frequency (E) variables. The Fonrier transform to frequency space can be 
reduced to an integral which can only be evaluated numerically. The results and the fits to it to a 
functional form for the imaginary part are given in Ref. (I22|) . We qnote this resnlt: 


x"ici,E,T) 


—Xo tanh 



q- 


1 

qF+^ 





( 3 ) 


(^) serves as an infra-red cut-off function due to deviations from criticality, where ke = 
replaces T as the infrared energy scale for critical fluctuations. Note that increases extremely 

)) from 0 on deviation from criticality. ^ 1 

at criticality when ke —t 0, so that the response is simply proportional to tanh(F^/2T) with an 
ultra-violet cut-off at Ec = For finite Kg, a fit to the numerical results gives. 


slowly, as an essential singularity (see Eq. 





—-for E/T <C 1; 

(l + ^keI2eT) 

1 (^1 + for EJT > E/T > 1. 


( 4 ) 


Kq = Eu {E/Ec) serves as an nltraviolet cutoff function Tu(0) = 1, J-u(cx)) = 0. In experiments, 
the ultra-violet cut-off scale may come from physics not in the effective low-energy model, as for 
example the Fermi-energy, if it is smaller than |rc|“^. We will be confronted with this situation 
below in the measurements in the heavy-fermion compound CeCng-xAua,. 

The most striking result in Eq. © is that x^^(q! -F) F) is a separable function of (E, T) and of q. 
Other remarkable results are that the dimensionless characteristic correlation wave-vector depends 
logarithmically on the characteristic correlation frequency. This means that the dynamical critical 
exponent is z = oo. Also, Eq. 0 for gives that the deviation of the characteristic correlation 
frequency with departure from criticality is determined by an essential singularity from its value 0 
at criticality. In other words, the cross-over from what is usually called the quantum-fluctuation 
regime, where the fluctuation scale is T to the truly quantum-regime, where the fluctuations are 
determined by {p — Pc)/Pc is extremely slow. 

The frequency and temperature dependence of these flnctnations at criticality is the same as that 
proposed phenomenologically [21] to give the marginal Eermi-liquid of fermions in the cuprates. The 
result for the spatial correlation length are quite different and remove the un-appetizing feature of 
the phenomenology that the spatial correlation length was Independent of deviation from criticality. 
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However, the single particle self-energy and the resistivity calculated due to scattering from the 
fluctuations given by Eq. © is also proportional to T in the quantum-critical regime, T/ke ^ 1 
(see supp. sec, in Ref. [TH]) and the entropy is proportional to TlnT. A microscopic origin m. 
quite different from AFM, has been found for criticality in the relevant range of cuprates, whose 
fluctuations also map to the 2d-XY model. 

The most stringent test of the theory is through the measurement of x(q) E, T), from which most 
other properties can be derived. These are very difficult measurements which also require large 
crystals. We compare here the results with the existing measurements in a Fe-based compound 
BaFei, 85 Coo.i 5 As 2 , and a heavy-Fermion compound CeCug-xAua;. These are the only measure¬ 
ments that we are aware of which are done in the frequency, momentum and temperature region 
of interest to test theories of AFM quantum-criticality. We will also recall some old results in a 
cuprate compound near the low doping where its AFM transition temperature appears —)• 0 but in 
which, due to the disorder, a spin-glass type phase sets in at a finite temperature. 


III. x"(q,A,T) IN BaFei.85Coo,i5As2 

BaFei. 85 Coo.i 5 As 2 has a putative Antiferromagnetic quantum critical point of the planar variety 
very close to the above composition, as we will show from the measurement of the spatial correlation 
length [25]. The fluctuations in the range of frequencies and temperature measured are about 
Q = (1/2,1/2,1) and equivalent points. The compound is strongly two-dimensional and the 
fluctuations are studied only in the {h, k)- plane. The superconducting transition temperature of 
25 K cuts off the low frequency part of the critical fluctuations and gives a peak in the frequency 
dependence at near twice the superconducting gap, about 10 meV, and half-width of peak is about 
4 meV. The q-dependence has been measured |26| at three different frequencies, 3, 9.5 and 16 
meV and at two different temperatures, one at 4 K, well below the superconducting transition 
temperature, and the other at 60 K. The frequency dependence at the critical vector for many 
frequencies from about 0.5 meV to 60 meV has been measured at these two temperatures and at 
280 K. We have to discard the 4 K data below about 15 meV, because of the superconductivity 
induced low energy features below about twice the superconducting gap, as further discussed in 
the caption to Fig. Q. A great virtue of the measurements is that the absolute intensity of the 
magnetic scattering has been measured. We therefore need to normalize all the data presented 
only once for all frequencies and momenta. 

We present the data for the q-dependence at the various indicated frequencies and temperatures 
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in Fig. Q. As shown, the distribution in q about the maximum fits a Lorentzian, with a width 
Kq ~ 0.04 ± 0.007 r.l.u., which is frequency and temperature independent to within the error bars, 
in the range (a factor of 5 in frequency and 15 in temperature) that it has been measured. This 
is consistent with the theoretical result that the q and the FI, T-dependence are separable. (The 
discarded data at 4 K for frequencies below twice the superconducting gap shows a 20% smaller 

Kq.) 



FIG. 1. x"(q, A, r) as functions of |q — Q| at various fixed temperatures and energies. Here, the q-scan 
is along (g,q, 1) and Q = (0.5,0.5,1). A background contribution to y" has been subtracted (assuming a 
linear dependence on q as given in Ref. [26] and fitted by the Lorentzian l/[(|q — QI/k^)^ + 1] (the curve 
shown has Kq=0.04 r.l.u.). 

Taking the measured large upper cut-off m of the fluctuations in this compound of about 200 
meV, and using the measured Kq k, 0.04rlu, Eq. Q gives ke ^ dAT, the lowest temperature mea¬ 
sured. Therefore, we should compare the frequency dependence with the form expected at critical¬ 
ity, i.e. oc tanh(E/2r). Fig. Q shows the absolute measurement of x"(Q,F;,r)/x"(Q,E/r»i) 
at the peak of its momentum-dependence at 60 K and 280 K for all values of energy measured. 
Also shown is the function tanh(£'/2T). In an inset we include the data at 4 K also. Above the 
superconductivity induced reduction followed by the bump, it also joins on to the data at 60 K 
and 280 K data in absolute value. This means that once we fix the magnitude xo of Eq. ([^, all 
the frequency and temperature dependent data is consistent with the form tanh(E/2T). 
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FIG. 2. x"(Q,-E,r) as functions of F;/(2T) with Q=(0.5,0.5,1) at r=60K and 280K for BaFei.gsCoo.isAss. 
The data is taken from Ref. [2S] . The inset shows the same data on a log-log scale, so as to also include the 
T = AK data (well below the superconducting transition temperature), showing that for Fi > 20 meV and 
up to 50 meV, the data is quantitatively is the same at the higher temperatures for E/2T ^ 1, as given by 
the theory. 

Two of the principal predictions of the theory are therefore shown to be obeyed. We urge 
measurements at various other compositions to test the dependence of the correlation lengths on 
deviations from criticality as well as more detailed E and T measurements. 


IV. x"(Q,F;, T) IN CeCu6_,Au, 

In CeCue-xAua:, Ising long range AFM order occurs for x > Xc = 0.1 at an incommensurate 
vector Qo ~ (0.6,0,0.3) [28j up to x = 0.3. The magnetic moments are aligned in the c-direction 
|28| . At the critical concentration Xc the AFM fluctuations are strongest at (0.6,0,0.3) but also 
at a wave vector (0.8,0,0) |29l I30j . which only develops short range order for higher x = 0.2 
|28j . Similar response E,T) for both Q vectors (and equivalent positions) were observed in 

the measured neutron scattering function for x ~ Xc along rods in reciprocal space [Ml El] from 
Qi = (0.8, 0,0) towards Qo = (0.6,0, ±0.3) as sketched in Fig.(|^). This is consistent with two- 
dimensional magnetic correlations in real space, perpendicular to the rods, along the b-axis and in 
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the (ac)-plane, as illustrated schematically in Fig. ([^) for one example. The neutron scattering 
data presented here was collected from different spectrometers [291 ED and calibrated using the 
incoherent nuclear scattering from the sample. 



FIG. 3. (a) The critical wave vectors Q of CeCus.gAuo.ij shown in reciprocal lattice. The experimental 
g-scan trajectories around Qi and Q 3 are also shown, (b) The plane of 2d fluctuations is shown in real 
space. Perpendicular to this plane, x"(q, F,r) is nearly independent of q, testifying the 2d nature of the 
fluctuation. Note that the plane of critical fluctuations is not a simple crystallographic plane. 

We first consider the q-dependence of the measured susceptibility [29] x" ^ with the fre¬ 
quency and the temperature fixed near two ordering wavevectors Qi=(0.8,0,0) and Q 3 =( 1 . 4 , 0 , 0 . 3 ). 
The q-scan directions are also shown in Fig. along 6 -axis for Qi, and along a trajectory in (ac)- 
plane for Q 3 . We normalize the results for different frequencies and temperature to their peak 
value at q = Q. Such normalized plots as a function of the deviation of q from Qi are shown in 
Fig. §r). A simple Lorentzian, as in Eq. ^ hts the data, with a width Kq ~ O.ISA which is 
independent of temperature or frequency in the range measured and within the error bars of the 
data. This independence is a test of the separability of the q and E-dependence of the fluctua¬ 
tions. In Fig. (|^), we show measurements near Q 3 at different temperatures m, which have 
much larger error bar (due to the low transfer energy). The best fit to the data shows the same 
width, independent of frequency and temperature as in Fig. &)• 

We now turn to the frequency and temperature dependence at q = Qi ^3 as well as Q 2 =( 1 . 2 , 0 , 0 ). 
From Kq Ri 0.13 A“^, we estimate that the deviation from criticality is very small for the sample 
measured, with -s/ {p — Pc) I Pc ~ 0.1. Taking E(. ~ 4iF, as will be justified below, we find from Eq. 
Q that ke ~ 4 X K, which is much smaller than the range of temperature measurements. 
Therefore, one is well justified in comparing with the theory for the dependence on frequency 
and temperature at criticality, Eq. ([^, with = 0) where it should simply be proportional to 
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FIG. 4. x"(q, i?, r) as functions of q for two g-scans around Qi=(0.8,0,0) (a) and Q3=(1.4, 0, 0.3) (b) at 
various fixed E and T for CeCu 5 . 9 Auo. 1 - A q-independent background contribution has been subtracted. 
The fitting curve is Lorentzian l/[{q — qcM /+ 1] with Kq=0.11 r.I.u. « 0.13 A“^ (considering b = S.lA) 
in (a) and Kq=0.13 A“^ in (b). 


tanh(i?/2T). However, unlike the case of the Fe compound discussed above, the E and T of the 
measurements go well across the Fermi-energy of about 4 K, estimated from the linear part of the 
specific heat [32]. The cut-off function TuiEjEM) due to the upper cutoff E^ can no longer be 
approximated as 1. The measured |29l El] ^'{Q„E,T), scaled by a constant xo [presumably the 
value of x^^(Q) E, T) at E/{2T) —)• 00 and E <C E^], is shown as a function of E/2T in Fig.([^). The 
data agrees reasonably with the function tanh(i?/2T) when E <C Ec, but systematically deviates 
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FIG. 5. x"(Q,i?,T) as functions of E/(2T) for various constant-i? or T scans for CeCus.gAuo.i at three 

ordering wave vectors Q = Qi, 2 , 3 - The label for different symbols shows the location of Q as well as the fixed 
value of E (or T), while T (or E) is varied. In (a), x"(Q: T) is scaled by a constant xo = 5.5fig{meV)~^ 
while in (b) it is scaled additionally by Eu{E/Ec) = l/y^l + {EjEcY with Ac = 4iF. 


when E > Ec ■ We now choose a cut-off function 

(k) ~ Vl + (E/E,)2 ■ 

We replot the same data, x"(Q>-^j^)/xo, divided by Eu{E/Ec), as a function of E/2T in Fig. 
&). With a value of Ec = 4K, the data, within its considerable error bars, is consistent with the 
scaling function tanh (E/2T). 

A complete test of the theory requires measurements varying x or pressure to vary the distance 
to criticality and thereby test the theory of the correlation length. It would also be desirable to 
have more measurements for the momentum-energy and temperature dependence for smaller E/T. 
CeCue-xAua; has a rather complicated magnetic structure. We urge neutron scattering results also 
on other heavy-fermions with simpler antiferromagnetic structure near their quantum-criticality. 

The q and E dependence of the data also has been fitted to an alternate phenomenological form 
earlier [31]. However, it does not give the observed linear in T resistivity. 


V. DISCUSSION 

The limited available data is consistent with the separability of the momentum and energy 
dependence of the critical fluctuations in two completely different experimental system, which 
share the feature that they both have itinerant AFM quantum critical points. The tanh(A/2T) 
dependence for energies smaller than the cut-off is also consistent with the data. As mentioned 
the linear in T-resistivity, and the TlnT entropy and thermopower are also properties of transport 
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due to coupling of fermions to such fluctuations. Since an external field provides a linear change to 
the energy for spin-systems, one expects the resistivity to be a function of {\H\/T). The resistance 
variation in a field at various temperatures has been found to be proportional to 

[1]. Magnetic fluctuations of the form of Eq. Q provide a constant contribution to the nuclear 
relaxation rate. This has also been observed [3] near quantum-criticality. 

Finally, we note that in cuprates for dopings near the quantum-critical region of the AFM but 
above a temperature, which appears to be determined by impurities, the resistivity is linear in T. 
In this region the spectral function of AFM-fluctuations, as determined by neutron scattering, does 
have E/T scaling and a temperature independent correlation length. See Refs. [33land[3^ 

There are many other materials showing quasi-2d fluctuations, and/or linear-T resistivity, which 
could also be tested by the criticality of 2d quantum-dissipative XY model. Heavy Fermions are 
notoriously hard to grow as large single crystals suitable for inelastic neutron scattering but the 
Ce compounds in the (115) family appear to be suitable for the purpose. Measurements have not 
been done on them in a suitable range of (q, E, T) to check the theory. The Fe superconductors 
appear also appear to be suitable for the purpose but systematic quantitative measurements of the 
kind needed are scarce. We urge such measurements. 
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